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Abstract 

The contributions to non-singlet evolution kernels P(z), for the DGLAP equation, and V(x,y), for 
the Brodsky-Lepage equation, are calculated for certain classes of diagrams including the renormalon 
chains. Closed expressions are obtained for the sums of contributions associated with these diagram 
classes. Calculations are performed in the [v? 3 ]e model and QCD in the MS scheme. The contribution 
of one of the classes of diagrams dominates for a large number of flavours Nf 1. For the latter case, 
a simple solution to the Brodsky-Lepage evolution equation is obtained. 
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1 Introduction 

Evolution equations play an important role in both inclusive H and exclusive [2] hard processes. They 
describe the dependence of parton distribution functions of DIS and of the parton wave functions on the 
renormalization parameter g 2 . The main ingredients of these equations are the perturbatively calcula¬ 
ble evolution kernels. The two-loop DGLAP evolution kernels P(z) were obtained in [p|], and a more 
complicated two-loop kernel V ( x , y ) of the Brodsky-Lepage (BL) evolution equation for meson wave 
functions was calculated in Si- A three-loop calculation of these kernels looks like being a tremendous 
problem. Recently, the results of very complicated calculations of the first few elements of the three-loop 
anomalous dimension of composite operators in DIS - 7^) (IV = 2,4, 6 , 8 ) have been presented in ||. 
Lately, these results were applied to improve the QCD analysis of DIS experimental data [j7[|. Despite 
this important technical and phenomenological progress one may feel insufficiency because of the absence 
of the whole kernel in this order or, even at least, some parts of it. Of course, the kernel admitting a 
physical probability interpretation is a more general and convenient object from the theoretical point of 
view. 

A method of estimating anomalous dimensions of composite operators in the limit of a large number 
of flavours, Nf, has been suggested by J. Gracey [j|, ||. It is based on conformal properties of the theory 
at the critical point g = g c corresponding to the non-trivial zero g c of the D-dimensional /3-function 
f3(g c ) = 0. The generating function has been constructed to obtain the leading larg e-Nf asymptotics of 
the anomalous dimension 7 ( n )(A) to any order n of perturbation theory (PT). 

In this letter, I propose a method of calculating certain classes of multiloop diagrams directly for the 
DGLAP kernel P(z) as well as for the BL kernel V(x,y) in the MS scheme rather than their moments 
as in Hi- These diagrams contain insertion of chains of one-loop self-energy parts (renormalon chains) 
into the one-loop diagrams (see Fig.l a,b) for kernels. The kernels P^ are obtained for any order n of 
PT by analysing of these “dressed by chain” diagrams. Then, the kernel P(z;g ) ( V ) for diagrams with 
the totally dressed propagator (i.e., in the “all loops” approximation) is calculated; this kernel appears to 

1 E-mail: miklis@thsunl.jinr.dubna.su 


1 



be a generating function to obtain the partial kernels P( n \. Moreover, this kernel P(z: g ) is the analytic 
function in variable of a coupling constant g. The insertions of the renormalon chains lead to converged 
PT series, and the usual renormalon singularities don’t appear in the renormgroup function calculations. 
This was found first for /3- and 7 -functions in the following development of this approach is presented 
in [11], It is to be stressed that the method of obtaining P or V does not depend on the nature of self¬ 
energy insertions and does not appeal to the value of parameters NjTr, Ca /2 or Cf (for QCD case) 
associated with different loops. Another distinctive feature is that the PT-improved evolution kernels 
are calculated in the direct and standard way by using rather elementary methods. This is possible due 
to a simple algebraic structure of the counterterms of the diagrams considered. In this letter, we mainly 
present technical results on resummation of the above mentioned counterterms within the framework of 
the [y ? 3 ]6 model and, in part, in QCD. The solution to the BL evolution equation for a large number of 
flavors Nf is also briefly discussed. 

To develop the diagrammatic analysis of multiloop evolution kernels, let us first consider a toy model 

describing ip*, ipi, the charged “quark” fields, and ip, the “gluon” one, in 6 space-time dimensions based 

N f 

on the Lagrangian Li nt = < 7 ^( 6 )- The number of “quark” flavors, Nf, will be considered as 


an arbitrary free parameter associated with the “quark” loops. This theory has much in common with 
QCD: it is renormalizable and its /1-function has a structure similar to that of (3qcd , 


/3(a) = -a 2 b 0 + 


u 0 = 


*1 

6 


a = 


(4vr) 


Moreover, the one-loop evolution kernels Po and Vo, following from the simplest triangular diagram in 
Figs.la and lb, turn out to be proportional to the corresponding expressions for the same diagrams in 
QCD 


( x x\ 

6(y > x)- + 0(y > x)- 

y yJ 


(i) 


where y = 1 — y, x = 1 — x,.... The similarity of the kernels in both the theories extends to the two-loop 
level [|l^]. The advantage of the [<£> 3 ]6 model is that the relevant Feynman integrals are simple, and it is 
much easier to study the structure of multiloop expressions. 

In general, calculation of the evolution kernels P(z ; a) [3| and V ( x, y, a) || is quite analogous to that 
of the anomalous dimensions. The major modification is the change ( kn) N —> S(z — kn) of the vertex 
factor corresponding to the composite operator, k being the relevant momentum related to the incoming 
“quark” line, see Figs. la,b. In these Figs., p ( yp or p) is the external momentum of the diagrams; n is 
a light-like vector (n 2 = 0) introduced to pick out the symmetric traceless composite operator; pn = 1. 


Detailed examples of similar calculations can be found in [12]. 


2 First triangular diagrams for DGLAP evolution kernel 


In this section, the method will be demonstrated by using well-understood [9, 10, O] quark-loop insertions 


as an example. Let us consider the one-loop triangular diagram T for kernels in Fig.lc with the chain 
of self-energy part insertions gi into the gluon line. Expressions for P(z) or V{x,y) in terms of the 
renormalization constant Zr for the diagram T in the MS scheme (see [|L4|] and ||) are given by 


Z r = 1 — KR\T), P = -ad a (z{, 1} ) = ad a (KiR'(T)) 


( 2 ) 


Here, R 1 is the incomplete BPHZ P-operation; e = (6 — D )/2 and D is the space-time dimension; K 
picks out poles in e, where K\ projects a simple pole, and Zd) is the coefficient of the simple pole in 
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the expansion of Z p. For g.;, being the one-loop insertions, expression (||) can be simplified following the 
definition of the -RhopCTafion as 


KiR'ir) = KiR'(G = Ki 




P(n ) = (n + l)aKi 


Ki)gi 


(3) 


where G denotes the “outer triangle” part of the diagram T without the chain of bubbles, and 

^ a,. n „2^ _ aN f (__ /V V X 


R9i(k 2 ) = (1 -K 1 )g i (k 2 ) = — 
7 v ,(e) = C(e)B(2 - e,2 - e). 


7*00 ^ - 7 *( 0 ) 


(4) 

(5) 


Here B{a,b ) is the Euler B-function, 7 ^ = 7 ^( 0 ) is the one-loop anomalous dimension of the gluon field 
(at Nf = 1), and the constant C'(e) = T(1 — e)r(l + e) reflects the concrete choice of MS scheme where 
every loop integral is multiplied by the scheme factor T(D/2 — l)(/r 2 /47r) £ . Note, the function 7 (e) in 
(j5j) will play an important role in our consideration. Substituting ([|) into (]|) and performing a direct 
calculation (see, e.g., [ 12 ]), one arrives at the expression 


P, 


(nj (^0 = ( n + l)aP 0 (z)(-aN flfip ) n K 1 


n-j 


C(e)( 1-e) A /7*(g)V r(l + V + l)e) (n\ (-) 

£ n+1 ~ ~P 0 \7*(0 )J r(l + je)Y(l + e) \jJ j + 1 


( 6 ) 


It is instructive to consider the properties of the first few kernels P$ (z). Their expressions are presented 
in Table 1 for k = 1,... 8. To obtain them, the FORM 2.0 program jl5| has been used f^. The lessons 
are : i) the perturbation theory for kernels (z) looks like being improved in comparison with the 

corresponding anomalous dimensions 7 ( n )(N). Indeed, the term in p/jj (z) leading in z is ~-he., 


>0 ; 


n! 


is not suppressed; 


(n) 

it is factorial suppressed in n, while for the corresponding contribution to 7 n (IV) ~ (n 
ii) one can see factorial suppression in n of all other logarithmic terms in P^Jj (z); iii) the ?r-bubble chain 

generates C( n ) (the Riemann zeta-function) in the non-logarithmic term in P^(z); this may be traced 
back to the expansion of the Euler B-function as their possible origin. These properties give hint about 
a possible resummation of the P^j (z)-series. 

OO 

Indeed, the closed expression for the sum of partial kernels P^\z]A) = a J2 P M (*) exist > 


n =0 


P (1 \z;A) =aJ2 p ln)( z ) = aP o( z ) z A (! - A ) ^ ^ 4 ) j ’ where A = aN Pv 


(7) 


and the kernel P^\z',A) is the generating function for P^jj (z). The proof of this assertion is presented 
in Appendix A. This result possesses several remarkable properties: 

• the Pd)(z; A) becomes the dominant part of the total perturbative kernel P(z) when 

JVj > 1. Below, the result for the kernel PW (z; A) will be completed by taking into account similar 
corrections to the “quark leg” (see Fig.ld) whose contribution is proportional to 5(1 — z). 

• the analytic properties of PP (z; A) in the variable A are completely determined by the one-loop 
“anomalous dimension” 7 ^(A) in D dimensions (see Eqs. (0 !))• The singularities of P^(z]A) 
correspond to zeros of the function B (2 — A, 2 — A) in 7 MA). Of course, 7 ^(A) provides the scheme 
dependence in PP(z; A) due to the factor C in the definition (||). 

2 I am greatly indebted to Dr. L. Avdeev who provided me with his brilliant FORM-based program for expansion in e, 
see ©■ Note that the contents of Table 1 is limited here only by place. 
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• the kernel P^\z]A) is an analytic function in variable A , except for singularities at the points 
A = 2 + k + 1/2, k = 0, 1,where the function has simple poles. The nearest singularity appears 
at A = 5/2, i.e., at aNf = 15, determining the range of convergence of the PT series. 

In a similar way, one can derive an expression for the sum of diagrams in Fig.l d connected with 
the anomalous dimension of the quark field 7 ^ (T) obtained in the “all quark-loop” approximation. 
Collecting the results of resummation in the main approximation in A, which correspond to the diagrams 
in Figs.lc and d, we arrive at the final expression 


P (1) (z;A)-5{l 



zz A (1 — A) — 5(1 — z) 


(3-A)(2-A)_ ' 


( 8 ) 


Integration produces the following expression for the anomalous dimension 7(IV, A) of the composite 
operator 


l(N,A) = J z N [p w (z-A)-6(l-z)^(A))dz = 

0 

r(4 - 2 A) 

a7 ^r (2 - A)r(i - A) 2 r(i + A) 

Formula (p) can be obtained by another method applied in [j|] to the QCD case (see Eq.(|27|) below). 
Note that the anomalous dimension 7 ( 1 , A ) corresponds to parton momentum. The equality 7 ( 1 , A) = 0, 
following from the r.h.s. of (9|), is expected; it signals that the parton momentum is conserved. 


l\Jy + l-A) F(2 -A) 


T(N + 3-A) r(4 — A) 


(9) 


3 Other triangular diagrams contributions 


Consider now the class of diagrams in Fig. 2. There are closed expressions for sums of partial kernels 


P^\z,B) = a Y,P™( z ) = aP 0 (z)(tz)- B ^A\ 


m =0 


, diag. in Fig.2a , (10) 

in Fig. 2a + MC, (11) 


P (2o) (z; B) = a£(2-5o, n )P ( t ) ) (z) = aPo(^)(^2(z)- B ^^-lj,diag. 

P (2b \z ; B) = a jr (m + 1 )P®)( Z ) = aP 0 (z) (l + B-^J > dia S- in Fi S' 2b > ( 12 ) 


■5(l-*b< 2) = aS(l - z) p A) , diagr. in Fig.2c ; 


(13) 


(2 \ 

the functions P(z ; B ) appear as generating functions for the corresponding partial kernels. Here (z) 
is the partial kernel with m insertions into one of the quark lines (see Fig.2a); 7 ^,(e) is the one-loop 
anomalous dimension of the quark field in P-dimension; for this model 7 ^(e) = 7 ¥ )(e), 7 ^ = 7 i/>( 0 ) = 
7 ^,( 0 ), B = MC denotes a mirror-conjugate diagram. The first Eq.(|lC)|) corresponds to the 

diagrams in Fig. 2a where the chain of quark self-energy parts is substituted only into the left quark line 
of the triangle. To prove it, one has to repeat the way similar to proof in Appendix A. Equation ( 0 ) 
corresponds to the sum of diagrams in Fig.2a and its MC diagrams. This result will be used to restore 
the corresponding kernel V(x, y ) in the next section. The analytic properties of the functions pl 2a ) (z; B) 
and p( 2a \z-,B) in the parameter B are the same as for the kernel P^(z;P); they are determined by 
the function 7 ^(B). Equation ( |l2| ) corresponds to substitutions of the chains into the both quark lines 
of the triangle in Fig.2b. At least, Eq.([bj|) corresponds to contributions to the anomalous dimension of 
the quark field from the diagram in Fig.2c. 
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The contribution p( 2 \z;B) will be suppressed in the parameter B in comparison with P^(z',A), 
if Nf P> 1, i.e., A S> B. Note, however, that the iV-moments of the kernels P^ 2a P - 7 ( 2 )(iV, P) 
decrease in N slowler than A) corresponding to the kernel P^\ Therefore, at sufficiently large 

N 7 ^(NjB) > 7 ^) (IV, A) for any fixed parameters A and B. 

The expression for the kernel P(n.,m)(^; A, B), corresponding to insertions of different one-loop parts 
in both gluon (re-bubble insertions) and quark (rre-self-energy part insertions) lines of the triangular 
diagram (see Fig. 2d), is obtained as well. This formula is similar to Eq.(| 6 j) for P( n \ , but looks more 
cumbersome and is not shown here. The partial kernels P(n,m)H A, B) can be obtained by using the 
FORM program, in general, for any given n and rre. For illustration, we demonstrate here the first 


nontrivial kernel aP) 


(id) 


P(i,i) — Pq(z)AB 


(1 

r , si 

2 20 1 

r . 51 

2 31 

r , 8i 

L. 51 

(2 

ln (^) + 3 

~ ~9~ + 2 

ln (^) + 3 

-Is- 2 

Hz) + 3 

Hz) + 3 


80 

+ -9 - 2C(2) 


in comparison with partial kernels of the same order in a, P(2,o) (= P( 2 )) following from the 


expression 


for p! 1 ) and P(o, 2 ) (= 3P^), from p( 2fc ) 


P(2,o) — Po(z)A 2 ( — 


ln (z) + 3 


20 


— 9 " ) > P (°’ 2 ) “ 3P H) B ( 


2 ! 


ln(z) + - 


To complete the section, we conclude that contributions to kernel P from any one-loop insertions into 
the lines of the triangular diagram are available now. 

4 Triangular diagrams for the Brodsky—Lepage evolution kernel 


Here, some partial results of the bubble resummation for the BL kernels V are presented. We obtain 
them as a “byproduct” of our previous results for the kernel DGLAP P. The results obtained may be 
used, in particular, to check the regular calculation of the kernel V in high orders of PT. 

Note, the diagrams for the BL kernel differ from the DGLAP diagrams only by the “exclusive” 
kinematics of the input momentum, compare the diagrams in Fig.la and lb. So, one can repeat the 
proof similar to that in Appendix A for this case. There is another far more elegant way - to use 
exact relations between V and P kernels for triangular diagrams which were established for any order 
of perturbation theory in [||. These relations work for both the [p^e model and QCD. I quote these 
propositions without proofs. 

Let the diagram in Fig.lc have a contribution to the DGLAP kernel in the form P(z) = p(z)+ 
5(1 — z) • C\ then its contribution to the BL kernel is 


V(x,y ) = C9(y > x) I 

J I 

where C = 1 + (x —>• x, y —> y). Substituting Eq. 
the expression for V ,yl \x,y) 


~dz + 5(y — x) ■ C, 


(14) 




6(y > x) 


for P^\z) into relation (14) we immediately derive 

(1-4) 


1 -A 


2 <% ^ (3 — A)(2 — A) 


(15) 


The later expression may independently be verified by other relations by reducing any V to P (5|, [D]] 
(V —> P reduction) 


VW (x, y- A) = C0(y > x)F(-\ A) - 5(y - x) ■ C(A) -»• P ( 1 ) (z; A) 

V 

where P^ (z- A) = 6(1 — z)z—F(z\ A) — <5(1 — z) ■ C, 


(16) 
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indeed, substituting ( |15| ) into (jl0) we return to the same Eq.(||) for P^\z). Moreover, the first term of 
the Taylor expansion of V ( ' 1 \x,y] A) in A coincides with the results of the two-loop calculation in [12|. 
The contribution should dominate for Nf S> 1 in the whole kernel V. and moreover, the function 
V^(x,y) possesses an important symmetry. Really, the function V(x,y;A) = V^\x,y; A) ■ ( yy) l ~ A 
is symmetric under the change x <->• y, V(x,y) = V(y,x). This symmetry allows us to obtain the 
eigenfunctions p n (x) of the equation 

1 

J V^\x,y,A)ijj n (y)dy = A)i/j n (x), (17) 

o 

i>n(y) ~ (yy) d ^ A ' > ~ 1 Cn^ A ' > 2 \y-y), here d^(A) = D a /2 - 1, D A = 6 - 2A, (18) 


d^p(A) is the effective dimension of the scalar field when the anomalous dimension is taking into account, 
and ci°\z) are the Gegenbauer polynomials of an order of a. This form of eigenfunctions {'i/j n (y)} as well 
as the x y symmetry of the function V(x, y) are the consequences of conformal symmetry conservation 
(see ]l 8 | for details) for the sum of diagrams under consideration Q Equation ([Tl]) is tightly connected 
with the BL evolution equation in a special case when the /3-function [3 = 0 . In this case, the BL equation 
is simplified, the variables are separated, and the evolution equation may be reduced to (jl7). The partial 
solution to the BL equation turns out to be proportional to ip n (y) in (|lg|). 

For the diagrams in Fig. 2a plus its MC diagrams, the relation between R( 2a ) and pO a ) is (see i) 


V {2a \x,y) = C9{y > x)^ 


-P^ a \x) + -P {2a \x) 
y y 


Ip( 2 a) f 

y \yJ 


(19) 


Substituting ([lT|) into ([H]) one arrives at the expression for V^ 2a \x,y) 


V^ 2a \x, y) =Ca9(y > x) 


~ x l ~ B x l ~ B 
—z 1 - 

y y 




( 20 ) 


which satisfies the same check as in the previous case. Of course, this formula does not cover the diagrams 
in Fig. 2b, where both the quark lines are dressed. 


5 The results for QCD evolution kernels 

The assertions about diagram contributions to P and V, similar to the presented above, are also valid for 
the QCD diagrams. Their proof does not contain essentially new elements, but looks rather cumbersome. 
The complete results for the QCD evolution kernels will be presented in a subsequent paper. Here, at 
first, the QCD results for triangular diagrams in Fig. 1 c,d in the Feynman gauge are discussed. Based 
on the proof in Appendix A in the QCD version, one can derive the result for the sum of diagrams, in 
Fig. lc, in QCD 


P^\z-,A) = a,C F 2z ■ (1 - P^jz- A - a,C F ■ <5(1 - z) ( ”)”) (1 ) A) - l) . (21) 

Oi 

Here, a s = —; n(e) is twice the contribution to a one-loop D-dimensional anomalous dimension of 
the gluon field and e = (4 — D)/2; n(0) is the contribution to the standard anomalous dimension; 
A = — a s n(0). It should be emphasized again that the form of Exp.(^) does not depend on the nature of 
self-energy insertions into the gluon line. One can use it for the resummation of both the quark ( rsj N f T R ) 

3 When this letter was considered by the Editorial Board, the work [|l9[ appeared, where a similar solution has been 
presentned for the case of QCD evolution kernel. 


6 













and gluon (~ Ca/ 2) loops. Substituting into the general formula the well-known expressions for 
11 (e) from the quark or gluon (the ghost loop is also added) loops 


n 9 (e) = -8N f T R B(D/2,D/2)C(e), 

n 9 (e) = ^ B (D/2-l,D/2-l)(j^^)C(e), 


one obtains the expression pjp 


(z; 5) for the quark-loop insertions 


if c) (*;<*) 


cl s C p2zz s 


(D q /2~l) 2 B(2,2) 

B(D q /2,D q /2)C(5) 


—a s C f<5(1 


*) ( _W_I_^ 

^\B(D q /2,D q /2)C(5)(D q /2-l) ) 


( 22 ) 

(23) 


(24) 


where D q = A — 25, 5 = — a s II 9 (0) = a s NfTp-, C(5) is the scheme factor, 

3 

and the expression P g lc \z',e ) for the gluon-loop insertions 


7 (lc) 0^) 


a s Cp2zz 6 


—a s Cp5( 1 


10 (D g /2 - 1 ) 2 {P g - 1) 

3 (3 D g - 2)B{D g /2 - 1, D g /2 - l)C'(e) 

_ z j™ _ (Dg ~ 1 ) _ 

l 3 (3 D g - 2)B(D g /2 - 1, D g /2 - l)C(e) (D g /2 - 1) 



(25) 


where D g = 4 — 2e, e = — a s II s (0) = — o, s Ca- ■ The expression (p5|) for P g lc \z; e) has no any singularity 

o 

in the parameter e < 0 due to the asymptotic freedom. 

Consider Eq.(|24j) for P q lc \z;5) in detail. By adding the contributions from the diagrams in Fig.Id 
to (^4|), one can find that the second term cancels in part and the final expression turns out to be the 
expected “plus form” (see, e.g. Q) 


P^ d \z-5) =a s C F 2 


\^-s (D q /2 - i) 2 n g (o) ~ 
n,(5) 


i.e 


, J 1 P^ (z; 5)dz = 7 j lc ’ d) (N = 0,5) = 0 (26) 


due to the current conservation. The analytic properties of P q lc ’ d \z] 5) in <5 are the same as for its scalar 
analogy pT) ( 2 ; A) (see ( 0 ) in sect.2); they are determined by the behavior of the function n g ((i) in 5, see 
Eq. (|22|) . The nearest singularity of P q ^(z;5) in 5 appears at a s Nf = 15/4. The moments of Eq. ( |2Gj ) 
agree with the corresponding part of the generating function in [j8|] 


7 ( q lc ’ d) (N,5) 


2 N{D q + N- 1) 
3(D q /2 + N)(D q /2-l + N) 


(DJ 2 ~ 1 ) 

B(D q /2,D q /2)C{5)D q 


(27) 


see the first term in Eq.(14) in [8|], and ref. (jl3| (note that our moments differ in sign from the definition 
of the anomalous dimension of composite operators used there). 

To complete the QCD calculations of the kernel pC) ; we need the contribution from the last diagram 
in Fig.le with the chain in the gluon line that is inserted into the composite operator. It can be obtained 
in a similar way, as in the previous QCD-calculations for Eq.(^) 


p( le \z] A) 


a s Cp 2 


f 2z l ~ A 11(0) \ 
\i-z n(A )J ’ 


(28) 
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and the expression automatically has the “plus form”, see [|5j for details. Collecting Eq.(26) and Eq.(pq), 
one easily arrives at the complete QCD expression for Pq l \z] S) 


P< j 1 \z-,S) = a s C F 2- 


zz " 5 


(1 - « 5) 2 + 


2z 


1-5 


1-Z 


n,(Q) 

iW 


(29) 


The contribution Pq 1 \z;8) is gauge invariant. The moments of Pq 1 \z',5 ) again agree with the complete 
generating function obtained in [||] (see Eq.(14) there). To obtain the contribution from both kinds of 
insertions into the gluon line one would take certain substitutions in (|^): n^e) + —> 11^(5) and 

e —> 5, here e = —a(II g (0) + II 9 (0)). Of course, for this case a final result must be gauge dependent. 


6 Conclusion 

The method for calculating certain classes of multiloop diagrams for the kernel P(z ) of the non-singlet 
DGLAP evolution equation is presented. These multiloop diagrams appear due to the insertion of chains 
of one-loop self-energy parts (renormalon chains) into the lines of the one-loop diagrams for the kernel. 
Closed expressions P( 1,2 \z,a 1 , 2 ) are found for sums of all the diagrams which belong to two of the 
diagram classes, see sect. 2-3. The corresponding proofs are based on a simple algebraic structure of 
the counterterms for the diagrams under consideration; this is outlined in Appendix A. Moreover, the 
kernels pW(z, a*) are generating functions for the partial kernels P^(z) which correspond to the n-bubble 

insertion. The contribution P^^(z;ai) (corresponding to one of the diagram classes) would dominate in 
the kernel for Nf 3> 1. The analytic properties of the function pW(z,aQ in the variable a, are briefly 
discussed. The expressions for partial kernels P(n, m ){z\ ai, 02 ) for the diagrams of a “mixed class”, in any 
order of perturbation theory can also be obtained by using the FORM program. 

The contributions V^' (x,y, a*) to the Brodsky-Lepage kernel are obtained in sect. 4 for the same 
classes of diagrams as a “byproduct” of the previous technique. When Nf 1, a special solution to the 
Brodsky-Lepage equation is derived. We emphasize that the method of calculating the evolution kernels 
pW or Rh) does not depend on the nature of self-energy insertions and does not appeal to the value of 
parameters NCa/2 or C F (for QCD case) associated with different loops. 

The method and results are exemplified with a simple [y? 3 ]6 model. It is clear that all the results 
obtained above in the framework of the scalar model have a wider meaning and can be applied to the 
QCD case. Here, in sect. 5, some QCD results are presented too; in particular, the kernel P) 1 (z; 5), that 
corresponds to the diagram dressed by the main quark-loop chain, is derived. The anomalous dimension 
~/q 1 (IV, 5) corresponding to this kernel agrees with the generating function obtained earlier by Gracey 
[|], The contribution Pg^(z;e) from the diagrams with the gluon-loop chain is derived in the same 
way. The same can be obtained for the mixed case when both kinds of bubbles (quark and gluon) are 
involved. The completed QCD case will be considered in detail in a subsequent paper. 
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A Appendix 

Here, formular is proved; all others resummation formulae in the text above can be obtained in a 
similar manner. 
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Proof. One could split the sum in square brackets, Eq.(pj), into two parts 


T(1 + (j + l)e) n\(-) n -i 


^ g(jj ttt - w*>+^ 


where S (n+1) (e) = 


n+1 


E ( p ( £ )> 


3=0 


i—i _ (n T l'' 

r(l + (J - l)e)r(l +e) l j , 


(A.l) 

(A.2) 


7 ( e ) 

Here, we use the notation: F(e) = , S( n+ 1 ) is the sum including “all combinatoric factors” of 

the l.h.s. of (A.l), and the last term in the r.h.s. of (A.l) may be associated with the j = —1 term in 
sum. The term S( n+ \j in the r.h.s. of (|A.1| ) does not contribute to the pole part, by virtue of Lemma: 

<S , (n+i)( £: ) < 0(e n+1 ). The expression for P^(z) can be derived by substituting decomposition ( A.l ) into 
(§) and using the above mentioned Lemma 


P$(z) = aP 0 (z)(-A) n k i 
= aPo(z)(A) n Ki 
>(i) 


-n+l ' 


z "(1 — e) ( C'(e)5( n+ i)(e) + 


(-n 

m)\ 


1 


T n+1 ' 


F(e) 


= aPo(z) 


(AY 


n\ 


J d n 


\ 

1 de n 

i m l 

/ 


(A.3) 


£— 0 


Now, let us take a sum of Pgj(z) over all n using Exp.(A.3); the differentials are naturally collected into 
the “shift argument operator” exp ( Ad e ) 


P {1) (z ; A) = aJ2 P (n) ( z ) = aF o( z ) ] ex P ( Ad e) 


,_ e (l -g)7y(0) 
7v>( £ ) 


(A.4) 


e=0> 


Finally, performing the shift of argument in the squar bracket in (|A.4|), we arrive at the expression 


p(1) (z;A) = a £ P$(z) = aP 0 (z)z~ A (l - A) , 

and the partial kernels (z) appear in the Taylor expansion of P^ A (z;A) in A by construction 
Lemma: 5( n+1 )(e) < 0{e n+1 ) for any analytic function F(e) at point e = 0. 

Let us consider the expansion of every element of the sum in 


(A.5) 


F{e) j 


r(l +je) 


F(e) T(1 + {j — l)e)T(l + e) 


(A.6) 


in powers of e up to e n . Any power e m of this expansion is accompanied by powers j\ where l < m, and 
a coefficient that does not depend on j. Therefore, the expansion of S( n+1 )(e) as a whole in the power 
series e m will generate coefficients of the powers which are composed only of the elements proportional 

toE/(-r +1 -^ n+lX 


3=0 


J 


|. All these elements are equal to 0 for m < n in virtue of the identity 


E/(-r w f n+lN 


= 0, if l < n + 1. 


(A.7) 


3=0 


Therefore, the power expansion of 5( n+ i)(e) can be started with the power higher than e n . Note, the 
similar trick with the identity (A.7) had been used in [10]. The proof implies an obvious generalization of 
the elements ( A.6 ) that can be constructed as superpositions of T-functions depending on j only through 
the arguments like 1 + je, 1 + (j — l)e,.... ■ 
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Table 1. The results of the P^(z\ A) calculations, the C( n ) is the Riemann zeta-function; note that 
C(2) and the Euler constant 7 e do not appear in this expansion. 
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Captions for Figures 


Fig-1 

The diagrams in Figs, la - Id are for the scalar model (and QCD); dashed line - for “gluons”, solid 
line - for “quarks”; the slash on the line denotes the delta function S(z — kn), see Fig.la; black circle 
denotes the sum of all one-loop chains; MC denotes the mirror-conjugate diagram. The diagram in 
fig. le is for QCD; solid line - for quarks, dotted line - for gluons. 


Fig.2 

The dashed circle in Figs.2 denotes the sum of chains of one-loop self-energy “quark” parts. 
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